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Abstract 



In the ac polarized Raman spectra of NaA^Os we have found anoma- 
lous phonon broadening, and an energy shift of the low-frequency mode as 
a function of the temperature. These effects are related to the breaking of 
translational symmetry, caused by electrical disorder that originates from the 
fluctuating nature of the V 4 ' 5+ valence state of vanadium. The structural cor- 
relation length, obtained from comparisons between the measured and calcu- 
lated Raman scattering spectra, diverges at T < 5 K, indicating the existence 
of the long-range charge order at very low temperatures, probably at T = 
K. 

PACS: 78.30.-j, 78.30.Ly, 64.60.Cn, 78.30.Hv 

The investigation of the physical phenomena associated to the low-dimensional mag- 
netic structures, like a spin-Peierls (SP) transition discovered in the inorganic compound 
GuGeOaB, are of great importance for better understanding of strong electron correlations. 
A very interesting interplay between spin and charge dynamics results in a phase transition 
discovered in NaX^OsiH. At higher temperatures (T > 34if) the magnetic susceptibility 
of NaX^Os is in excellent agreement with the Bonner-Fisher curve for the one-dimensional 
Heisenberg antiferro magnet. At low temperatures, however, the susceptibility decreases 
rapidly to zercll suggesting a SP transition at T=34 K. On the other hand, the temperature 
dependent nuclear magnetic resonance (NMR) spectra^ showed a change of the vanadium 
valence across the phase transition, from uniform V 4,5+ to two different V 4+ and V 5+ states. 
These measurement sBq gave direct evidence for the charge ordering (CO) phase-transition 
scenario in NaX^Qs. A structural analysis,! and the polarized Raman and infrared (IR) 
spectra of NaX^O^, are also consistent with the existence of uniform vanadium valence in 
the high-temperature phase. However, despite extensive experimental and theoretical effort, 
no consistent picture has yet emerged for the low-temperature phase@~0. None of the mod- 
els proposed to date, have explained fully the two intimately connected issues: the nature 
of the phase transition and the low-temperature ground state. Raman spectroscopy is a 
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powerful method for the study of dynamics of solids and it can be used to address these 
issues. Again, despite a considerable amount of data, the Raman spectra of NaX^Os are still 
not completely understoodlB. This holds for the vibrational modes as well as for possible 
magnetic excitations. 

In this report, we present the study of the vibrational modes in the Raman spectra of 
NaV2C>5 and analysis of the effects associated with non conservation of the light-scattering 
wave vector selection rule k=0. We argue that anomalous phonon broadening and frequency 
shifts, that are observed in the ac spectra as a function of the temperature, are caused 
by strong electrical disorder due to the fluctuating nature of the V 4 5+ valence state of 
vanadium. The temperature dependence of the structural correlation length is obtained from 
the comparison between measured and calculated Raman scattering spectra, and used to 
characterize the phase transition in NaX^Os. This analysis shows, that the phase transition 
at T c = 34K corresponds to the onset of the short-range electron correlations, with a true 
long-range order (CO) at T ~0. Thus, the charge order is not static below Tc, in contrast 
to the conclusions of Lohmann et afi 

Polarized Raman scattering experiments were performed on NaX^Os single crystals (size 
~ 1 x 3 x 1 mm 3 along a, b, and c respectively) prepared as described in Refs. 2 and 3. As 
an excitation source we used the 514.5 nm laser line from an Ar + ion laser. The beam, with 
an average power of 5 mW, was focused (spot diameter ~ 80/mi) on the (001) and (101) 
surfaces of the crystals. The spectra were measured in a backscattering geometry using a 
dilor triple monochromator equipped with a LN2 cooled CCD camera. NaX^Os crystallizes 
above T c in the orthorombic centrosymmetric space group Pmmn (D^), with two molecules 
in the unit cell of a size: a=1.1318 nm, b=0.3611nm and c=0.4797 nm. Each vanadium 
atom, in the average valence state 4.5+, is surrounded by five oxygen atoms forming VO5 
pyramids. These pyramids are mutually connected via common edges and corners to form 
layers in the ab plane, see Fig. 1. The Na atoms are situated between these layers as 
intercalants. The structure of NaX^Os can also be described as an array of parallel ladders 
(running along the b direction) coupled in a trellis lattice (each rung is made of a V-O-V 
bond). 

The vibrational properties of NaVoOs are studied in great details by measuring the Ra- 
man spectra only from (001) surfacei'0. In our previous paper0 we observed an interesting 
low-frequency structure in the ac scattering geometry which, as we will show below, ex- 
hibits an unusual temperature dependence. The polarized Raman scattering spectra from 
(001) and (101) planes, at temperatures above and below phase-transition temperature, are 
presented in Fig. 2. The b axis of the crystals is set to be parallel to the laboratory H 
horizontal) axis. Thus the VV (V equals vertical) polarized scattering geometry from the 
(101) plane gives the mixture of aa, cc, and ac contributions (upper spectra in Fig. 2). The 
VV polarized configuration from the (001) plane gives only aa contribution (lower spectra 
in Fig. 2). The corresponding HH spectra (gives the bb contribution for both planes) are 
found to be identical. In this way we were able to determine the ac contribution at very low 
frequencies, otherwise impossible to get. For example, the low quality of the (010) surface, 
prevent us from direct measurements of the low frequency scattering in the ac geometry. 

In addition to the phonon modes, and a continuum centered at 650 cm -1 which belongs 
to the aa channel, we found three modes in the VV spectra obtained from the (101) planes 
(Fig. 2), that represent the ac contribution. Their detailed temperature dependencies are 



2 



shown in inset of Fig. 2. The lowest frequency ac mode is centered at about 107 cm -1 . 
The mode is asymmetric with a frequency cutoff around 120 cm -1 (Fig. 2 and Fig. 3b). As 
the temperature is decreased below T c , the mode becomes more symmetric, and hardens by 
amount of 10 cm -1 [Fig. 2 and Fig. 3b]. The integrated intensity of this feature increases 
rapidly by the increasing temperature below T c , and becomes constant for T > T c (Fig. 4). 
The temperature behavior of the frequency and the integrated intensity is similar to what is 
usually observed for a two-magnon modeE3. However, the two-magnon origin of this mode 
can be ruled out since its energy (117 cm -1 ) below the phase transition temperaturesE! is 
much smaller then 2A S ~ 160 cm -1 . The additional argument for excluding the two-magnon 
scattering process comes from observation of the mode at temperatures as high as 10 x T~. 
At these temperatures the two-magnon mode should not be visible in the Raman spectraE3. 
A one-magnon scattering process, and other magnetic-related scattering mechanisms, can 
be also eliminated since we did not find any change in the spectra in the magnetic fields up 
to 12 T. 

Our analysis shows that this structure might be related to the low frequency phonon 
density of states, due to the breakdown of the k=0 conservation rule in the light scatter- 
ing process. In NaX^Os the V 4,5+ valence state of vanadium causes the random nature of 
the coupling between the atomic displacements and the fluctuations of the dielectric sus- 
ceptibility. Such randomness actually appears because of irregular atomic bonding. The 
visualization of randomness in NaVoOs is schematically presented in Fig. 1 and it can be 
referred as an "electrical disorder" ,E3. Imagine that we "freeze" the electrical configuration 
of the charges in the rungs of the ladders, at certain time in the high-temperature phase. 
We find random electron configurations among rungs, Fig. lb. The average value for all rung 
configurations gives the V 4 5+ in which the electron is shared by the V atom at the each end 
of a rung. Below the phase transition temperature the charges start to order, (V 4,5±a , a^O) 
Fig. lc, reaching complete "zig-zag"0 order at T=0 (a = 0.5), Fig. Id. The "zig- zag" 
phase is presented as a real low-temperature geometry, even though there are some other 
proposed CO configurations!!!. It will be evident from our results that we were not able to 
discriminate among the various possible CO patterns, and Fig.l must be regarded only as 
illustrative of the electrical disorder. However, the "zig-zag" charge order is consistent with 
observed phonons in the Raman and IR spectra of NaV2Os0. If so, even for the perfect 
plane-wave phonon, for example an acoustic or optical phonon with a finite wave vector k, 
disorder of the atomic coupling allows inelastic scattering of light from this modeil. Then, 
the light scattering is expected to be proportional to the phonon density of states properly 
weighted by the coupling constant which is in fact a function of u. 

The light-scattering cross section is proportional to the Fourier transform of the correla- 
tion function of the polarizability fluctuations: 

h,M ~ / / dtd( ri - r 2 ) e -t-*(n-r 2 ) < 5x *(r u t)6 Xj (v 2 ,0) >, (1) 

where k = k/ — kg and uo = uj — u>s are wavevector and frequency of phonon(s) that 
participate in the light scattering process. The < ... > denotes the equilibrium expectation 
value. 

For the ideal crystals, both energy and wavevector conservation rules are fulfilled, and 
the first-order Raman intensity of phonons is proportional to two delta functions, 5(cjj — 
u) s -u) -5(k-q). 
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In disordered crystals the k=0 selection rule is broken due to lack of the translational 
invariance, and the Stokes part of the Raman scattering intensity is + n(u)] ~ 

J2j Cj{uj)gj{oj). The gj(u) is a density of states for the band j, and the n(u) is Bose distri- 
bution functional. Cj (o;ldescribes the coupling of the light and the vibrational mode u. In 
a continuum descriptions the fluctuations of the susceptibility result from the elastic strain 
field ej(u,r) = kej(uj)exp(ik.r) of the phonon (plane wave), which is modulated by static 
fluctuations 5pik(r) of the elasto-optical constants: <5xi( r ) — — e 2 / (47r) [p^- + 5pij(r)]ej(r). In 
this way, the coupling constant C(u) ~ k 2 J drexp(ihr)F(r) is expressed as a correlation 
function F(r) of the fluctuations bp of the elastoptical constants, which characterize the 
electrical disorder. The form of the correlations is usually postulated to be either exponen- 
tial damping exp(—r/l c ); or Gaussian damping exp(—r 2 /l 2 ), where l c is correlation length. 
In the case of NaX^Os we assume Gaussian damping to describe the electrical correlations. 
The correlation length may be defined as a length over which electrons in neighboring rungs 
"see" each other. In fact, this is just a positional (structural) correlation length. One sim- 
plified picture, where the intersite Coulomb interactions are "switched on" at T=34 K has 
been previously suggested^. 

Thus, the coupling constant is C(u) ~ k 2 e~ k lc ^. The to dependence of C(u) comes from 
the dispersion relation between the frequency and the wavevector [for example, if u = ck 
then C{oj) = uj / c) 2 exp[—{uj / c) 2 l 2 /A\. The same type of the correlation function and the 
coupling constant have been obtained by Martin and Brenigll in their analysis of Brillouin 
scattering in the amorphous solids. 

Finally, the normalized Raman intensity is: 

/(^) IT ^y - E k 2 e~^l%u - w (k)], (2) 

where a>(k) is a phonon dispersion. If we confine our analysis to the energy range of the 
acoustic or/and low frequency optical phonons (u < 150 cm -1 ) the Raman spectrum in the 
high-T phase is influenced by two contributions: acoustic or optic phonon density of state 
contribution, Y.k^ ~ w (k)] and the coupling function C{k). Since dispersion curves of 
the phonons have not yet been measured, we are forced to make assumption about which 
phonons are involved in the light scattering process. There are two possibilities: 

(a) Broad feature corresponds to the acoustic phonon with a zone-boundary energy of 
117 cm -1 . Bellow T c this mode is introduced to k=0 by zone folding effect. The x-ray 
diffraction experiments showed the existence of superlattice reflections below T c with a 
lattice modulation vector q=(l/2, 1/2, 1/4)0. In this case the strong anomaly of the elastic 
constants is expected at T c and indeed observed by Schwenk et al£3. 

(b) The mode corresponds to Iow-cj optical phonon, also allowed by the symmetry of 
the low-T phase, with energy that decreases as a function of the wavevector. According 
to the lattice dynamical calculations^ a good candidate for that phonon could be the low- 
uj B 2g phonon (active in ac polarized geometry), with Na vibrations mainly along the a 
axis. By examining the Raman spectra it is difficult to conclude which one of these two 
assumptions is more appropriate, because of the strong quasi elastic background at low 
frequencies. However, this choice does not critically influence our conclusions, and for the 
sake of simplicity we assume the cos k form of the phonon dispersion, uo = ujQCOsk/2, ujq = 117 
cm' 1 , k G [0, n/a]. 
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The calculated Raman spectra are presented in Fig. 3a and compared with measured 
ones, Fig 3b. Both the shift and the broadening of the mode are in good agreement with the 
experiment. The spectra are obtained by varying just one parameter l c , evaluating equation 
2 in one dimension, and by taking the values of k and l c in appropriate units of a-lattice 
constant. By increasing temperature (from T=0), disorder is introduced and the contribu- 
tion of the C(k) in the spectra becomes important. The increase of the degree of disorder 
is produced by the increase of the k ^ contributions, directed with C(k). Therefore, 
the broadening of the mode and its energy shift towards lower energies are produced by 
decreasing l c which is a measure of the degree of disorder; complete disorder is characterized 
with l c ~ a-interatomic distance and long range-order with l c — oo. The long-range order 
solution, l c = oo and C(k) = of equation 2, gives the vanishing of the Raman intensity. 
This is not unphysical. It is just telling us that one has to analyze the Raman spectra at 
T=0 using equation 1 instead of equation 2. 

The peak position of the mode as a function of the temperature is shown in Fig. 4a. The 
circles (full lines) represent the experiment (theory). Since in our calculation of the Raman 
spectra the temperature does not enter as a parameter, the temperature dependence of the 
peak maximum is included only through the temperature dependence of the correlation 
length. It is generally expected that correlation length changes with temperature, and 
this dependence can be obtained by examining the phonon frequency change close to the 
phase transition-temperature. The best agreement between the measured and the calculated 
spectra is obtained by assuming a quadratic relation between the correlation length and 
the temperature, l/l c ~ T 2 . Therefore, the peak position can be used as a measure of 
disorder. Accordingly, a similar temperature dependence is expected for the order parameter 
a. Furthermore, we present the comparison between the temperature dependence of the 
measured and calculated integrated intensities, Fig. 4.b. Here, we also find very good 
agreement between theory and experiment using the same quadratic relation between l/l c 
and T, and by introducing a non zero intensity offset at T=0 K (I(T = 0)/ I(T = 100) = 0.3). 

The temperature dependence of the correlation length is presented in inset of Fig. 4. 
The particular value of lc that is used to calculate the Raman spectra, does not have deeper 
physical meaning because of an arbitrary factor in the exponent of equation 2. But its tem- 
perature dependence does. The correlation length has a approximately constant value above 
T c and increases below T c . This signals that the critical temperature T c = 34 K represents 
the onset of short-range electron correlations. As the temperature is lowered below T c , the 
correlation length rapidly increases indicating the existence of a possible singularity at tem- 
peratures below T=5 K. In this case the divergence of the correlation length corresponds to 
the appearance of the true long-range charge order at very low temperatures (probably at 
T=0). The short range electron correlations could correspond to intersite Coulomb inter- 
action effectsii, which become important below 34 K, but the electron correlations should 
also persist in some form at temperatures above 34 K. Following the same arguments, the 
change of the spectral shape above 100 K (it is also found that some IR spectral changes 
occurs around 100 K,@) could be a consequence of an additional phase transition, magnetic 
in origin for example. 

In conclusion, we have presented the evidence for the existence of translational symmetry 
breaking effects in the Raman spectra of NaX^Os. Non conservation of the light-scattering 
wave vector selection rule, k « 0, is caused by strong electrical disorder due to the fluctuating 
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nature of the valence of vanadium. The temperature dependence of the structural correlation 
length has been obtained from the comparison between the measured and the calculated 
Raman scattering spectra. This suggests that the phase transition at T c = 34K, represents 
the onset of the short range electron correlations, with true long-range order developing at 
T ~0. 
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FIGURES 



FIG. 1. a) Schematic representation of the NaV20s crystal structure in the (001) plane. The full 
circles represent oxigens, the small (large) circles represent vanadium (sodium) atoms. Schematic 
representation of the electrical disorder (order) at b) T > T c , c) T < T c , d) T=0 in the vana- 
dium-oxygen ladders of NaV20s. Each rung in the ladders represents the V-O-V bond with one 
electron (black dot). 

FIG. 2. The polarized Raman scattering spectra of NaV20s at 300 K (thin line) and at 10 K 
(thick line), measured from (001) surface (lower spectra) and (101) surface (upper spectra). Arrows 
represent modes associated with the ac polarized contribution, see text. Inset: The temperature 
dependent (aa+cc+ac) Raman scattering spectra. 

FIG. 3. The a) calculated and b) measured Raman scattering spectra in 50 to 130 cm -1 
frequency range at various temperatures. 

FIG. 4. The a) frequency and b) normalized integrated intensity of the 117 cm -1 mode as a 
function of the temperature. The circles (full lines) represent the measurements (theory). The 
dashed line is guide to an eye. Inset: The temperature dependence of the correlation length. 
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M.J.Konstantinovic et al. Fig. 1 
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